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Abstract 

Background: Excision repair cross-complimentary group 1 (ERCCl) is an essential component of tlie nucleotide excision 
repair system that is responsible for repairing damaged DNA. Functional genetic variations in the ERCCl gene may alter DNA 
repair capacity and modulate cancer risk. The putative roles of ERCCl gene polymorphisms in lung cancer susceptibility 
have been widely investigated. However, the results remain controversial. 

Objectives: An updated meta-analysis was conducted to explore whether lung cancer risk could be attributed to the 
following ERCCl polymorphisms: rsll615 (T>C), rs3212986 (OA), rs3212961 (A>q, rs3212948 (OC), rs2298881 (OA). 

/Wef/jocfe; Several major databases (MEDLINE, EMBASE and Scopus) and the Chinese Biomedical database were searched for 
eligible studies. Crude odds ratios (ORs) with 95% confidence intervals (CIs) were used to measure the strength of 
associations. 

Results: S\yiXeen studies with 10,106 cases and 13,238 controls were included in this meta-analysis. Pooled ORs from 11 
eligible studies (8,215 cases vs. 1 1,402 controls) suggested a significant association of ERCCl rsl 1615 with increased risk for 
lung cancer (homozygous: CC versus TT, OR=1.24, 95% CI: 1.04-1.48, P = 0.02). However, such an association was 
disproportionately driven by a single study. Removal of that study led to null association. Moreover, initial analyses 
suggested that ERCCl rsl 1615 exerts a more profound effect on the susceptibility of non-smokers to lung cancer than that 
of smokers. Moreover, no statistically significant association was found between remaining ERCCl polymorphisms of 
interest and lung cancer risk, except for rs3212948 variation (heterozygous: CG vs.GG, OR = 0.78, 95% CI: 0.67-0.90, P = 0.001; 
dominant: CG/CC vs.GG, OR = 0.79, 95% CI: 0.69-0.91, P = 0.001). 

Conclusion: OyeraW, this meta-analysis suggests that ERCCl rs3212948 G>C, but not others, is a lung cancer risk-associated 
polymorphism. Carefully designed studies with large sample size involving different ethnicity, smoking status, and cancer 
types are needed to validate these findings. 
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Introduction 

It has become clear that the susceptibility to disease varies from 
one individual to another. Some heritable characteristics that 
shape the effects of environmental exposure may contribute to the 
variable susceptibilities among people. Smoking highly increases 
risk of lung cancer by up to 20 folds since many carcinogens in 
cigarette smoke can be converted into reactive metabolites in 
human body. These reactive products (e.g., diol epoxide deriva- 



tives of polycyclic aromatic hydrocarbons) can potentially damage 
cellular DNA and cause the formation of DNA adducts via 
covalent binding or oxidation. The resulting adducts are 
carcinogenic and may block the transcription of critical genes or 
result in mutations at hot spots [1] . Nevertheless, only a proportion 
of cigarette smokers and victims of second-hand smoke develop 
lung cancer in their lifetime, which may be partially attributed to 
the fact that DNA repair systems (e.g., nucleotide excision repair) 
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can effectively remove DNA lesion and restore genomic integrity. 
It has been suggested that DNA repair capacity is essential in 
protecting people from cigarette smoke-related carcinogenesis 
[2,3]. The nucleotide excision repair (NER) system is responsible 
for the repair of various DNA lesions, such as bulky adducts, cross 
links, oxidative DNA damage, thymidine dimers and alkylating 
damage. Genetic variations in DNA repair genes may alter DNA 
repair capacity and modulate cancer risk in the host. Excision 
repair cross-complimentary group 1 (ERCCl) is a critical protein 
in the NER pathway. Typically, ERCC 1 joins XPF endonuclease 
(also known as ERCC4) to form heterodimeric endonuclease 
piPF-ERCCl), which excises the 5' end of DNA to the damaged 
site. XPF-ERCGl complex also participates in homologous 
recombination and the repair of inter-strand crosslinks. Thus, 
functional polymorphisms in the ERCCl gene that compromise 
DNA repair capacity may be a potential risk factor for tobacco- 
induced cancers. Therefore, potential associations between ERCCl 
gene polymorphisms and cancer risk have evoked great interest. 
Due to lack of non-synonymous single nucleotide polymorphism 
(SNP) in the coding region of the ERCCl gene, most studies have 
focused on the rs3212986 (3'UTR C8092A) and synonymous 
rsll6I,T (exon 4 T190()7C) polymorphisms, which are believed to 
influence transcript stability and mRNA levels [4,,5], respectively. 
Numerous studies on such topics have been conducted over the 
past decades, but the results remain controversial. Several meta- 
analyses have also yielded conflicting conclusions [6-8]. Lately, 
more case-control studies regarding these topics have emerged. 
Apart from those two common ERCCl variants, associations of 
tiiree other ERCCl polymorphisms (rs3212961 (17677A>C), 
rs3212948 G>C, and rs2298881 OA) with lung cancer risk 
have also gained increasing attentions. Therefore, we conducted 
this updated meta-analysis to reassess the associations of the first 
two common polymorphisms in the ERCCl gene and lung cancer 
risk, and explore the influence of the other three polymorphisms 
on predisposition to lung cancer. 

Materials and Methods 

Publication search 

A systematic literature search throughout the MEDLINE, 
EMBASE and Scopus databases was performed with the use of the 
following search terms: "ERCCl or excision repair cross-compli- 
mentary group 1", "DNA repair", "polymorphism or variant or 
variation" and "lung cancer or tumor or carcinoma or neoplasm". 
Publications written in Chinese were also searched from the 
Ghinese Biomedical (GBM) database (http:/ /cbmwww.imicams. 
ae.cn/cbmbin) (1978—) to increase the coverage of our current 
study using the combination of terminologies: "ERCCl" , "poly- 
morphism" and "lung c:ancer" in Chinese. In addition, references 
of the retrieved research or review articles on this topic were 
manually reviewed to identify extra eligible studies. Publication 
search was initiated on October 12, 2013, and last search was 
performed on November 15, 2013. 

Selection of eligible studies 

Selected studies must meet the following inclusion criteria: 1) 
studies investigating any type or mixed t^pes of lung cancer (small 
cell lung cancer or non- small cell lung cancer (NSCLC): 
adenocarcinoma, squamous carcinoma, and large-cell carcinoma, 
etc.), regardless of smoking status; 2) original research studies 
written in English or Chinese, 3) case-control, nested case-control, 
or cohort study, and 4) adequate information to calculate odds 
ratio (OR) and 95% confidence interval (CI). Exclusion criteria 
were: 1) duplicate data, 2) abstract, case report, comment, review, 



and editorial, 3) lack of sufficient genotyping data, and 4) studies 
involving subjects with family history or cancer-prone disposition. 
In the case of studies with overlapping subjects, the latest and/ or 

largest study was chosen. Studies in which the genotype 
frequencies of the control group departed from Hardy-Weinberg 
equilibrium (HWE) was excluded from the final analysis, unless 
there was further evidence validating HWE from another ERCCl 
polymorphisms. 

Data Extraction 

Information extraction was conducted individually by two 
investigators (Zhu JH and Hua RX) from each eUgible study. The 
main information comprises first author, year of publication, 
countr)' of origin, ethnicity, cancer type, source of controls (i.e., 
population-based or hospital-based), genotyping method, number 
of cases and controls, genotype counts of five ERCCl polymor- 
phisms for cases and controls, and main findings. In the case that 
studies included subjects of different ethnic groups, data were 
obtained separately for each ethnic group and labeled as 
Caucasian or Asian. 

ERCCl gene expression analysis based on ERCCl variant 
genotypes 

The biological plausibility of our findings was investigated by 
correlating respective ERCCl polymorphism genotypes and 
corresponding ERCCl mRNA expression levels in 270 lympho- 
blastoid cell lines. Genotype data of ERCCl polymorphisms and 
ERCCl mRNA expression information were retrieved from 
HapMap website (http://www.hapmap.org) and SNPexp online 
tool (http://app3.titan.uio.no/biotools/help.php7app = snpexp), 
respectively. These resources facilitate researchers to determine 
the correlation of HapMap genotypes in a genomic region of 
interest and gene expression levels [9] . The international HapMap 
phase (II+III) release #28 data set contain genot\pe data of 3.96 
million polymorphisms for 270 individuals from four worldwide 
populations [CEU: 90 Utah residents with ancestry from northern 
and western Europe; CHB: 45 unrelated Han Chinese in Beijing; 
JPT: 45 unrelated Japanese in Tokyo; YRI: 90 Yoruba in Ibadan, 
Nigeria] [10]. The mRNA expression data were obtained from the 
same 270 individuals (GENe Expression VARiation, http://www. 
Sanger, ac.uk/ resources/ software/ genevar/). 

Statistical Methods 

Crude odds ratio (OR) with 95% confidence interval (CI) were 
used to access the association between each of the five ERCCl 
polymorphisms and lung cancer risk. Letters V and W represented 
variant and wild alleles of each polymorphism. Different genetic 
models were adopted during risk estimation: homozygous (W vs. 
WW), heterozygous (WV vs. WW), recessive (W vs. WV/WW), 
and dominant (W/WV vs. WW) model. A chi-square-based Q; 
test was used to check heterogeneity assumption. A P-value of > 
0.10 for the Q^test suggested the absence of heterogeneity across 
the studies, and then the fixed-effects model (the Mantel-Haenszel 
method) [1 1] would be adopted; otherwise, the random-effects 
model (the DerSimonian and Laird method) [12] would be 
performed. The Z-test was used to determine the significance of 
the pooled OR. The funnel plot was created to test publication 
bias. Briefly, the standard error of log (OR) of each inv(;stigati()n 
was plotted against its log (OR), and the asymmetry of funnel plot 
was assessed by the method of Egger's linear regression test [13]. 
Furthermore, sensitivity analysis was used to determine stability of 
the result, i.e., an individual study was excluded at a time, and 
then risk estimates were recalculated to examine the effect of single 
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study on the pooled ORs. Hardy-Weinberg equilibrium (HWE) in 
the control group for each ERCCl polymorphism was checked by 
the Pearson's goodness-of-fit chi-square test. Student's t test and 
analysis of variance test were used to evaluate the differences in the 
relative mRNA expression levels among different genotype groups. 
All analyses were performed by using STATA version 1 1.0 (Stata 
Corporation, College Station, TX) and SAS version 9.1 (SAS 
Institute, Cary, NC). All statistical analyses were two-sided, and 
P<0.05 was considered significant. 

Results 

Characteristics of studies 

Literature search initially produced 129 articles. After screening 
and reviewing articles for study eligibility, 113 articles were 
excluded due to review articles or non case-control studies, 
overlapped participants, failure to report genotype distribution 
data, and other reasons destxibed previously [14] (Figure 1). For 
example, among the 15 studies concerning ERCCl rsll615 
polymorphism [4,15-28], 5 were reported by Yin JY et 
a/. [18,25-28], in which samples were collected in the same 
institute, Liaoning Cancer Hospital. One of them [25] was not a 
case-control study and thus was excluded. In order to avoid 
repetition in the sampling, we chose one [18] from the- rest of 4 
studies [18,26-28], which has the largest sample siz(- and also was 
published most lately. Ultimately, 16 studies [4,5,15-24,28-31] 
with 10,106 cases and 13,238 controls were chosen for this meta- 
analysis (Figure 1). Of the resulting 16 studies, 9 were conducted in 
Chinese population, and 7 in Caucasian population. All studies 
report(;d the association between at least one ERCCl polymor- 
phism of interest and lung cancer risk. Studies, inspecting the 
associations between multiple ERCCl genetic variations and lung 
cancer risk were divided into several sub-studies, each of which 
encompassed the analysis of a single polymorphism. Genotypingin 
these 16 studies was performed via different methods, details of 
which were listed in Table 1. Overall, the current meta-analysis 
contained 11 studies for rsll615 [4,15-24], 6 for rs32 12986 
[4,15,16,20,22,31], 4 for rs32 12961 [4,16,28,30], 3 for rs3212948 
[5,29,30], and 4 for rs2298881 [20,22,28,29] polymorphisms, 
respectively. Most of these studies concerned all t^pes of lung 
cancer among both smokers and non-smokers, except for two. 
One study conducted by Zienolddiny et al. [4] focused on NSCLC 
in smokers and former smokers, while the other by Yin et al. [17] 
investigated only adenocarcinoma in non-smokers. 

ERCCl rs11615 (T>C) polymorphism 

In the meta-analysis of 11 studies with 8,215 cases and 11,402 
controls, the pooled OR for the association between ERCCl 
rslI615 polymorphism and lung cancer risk was statistically 
significant (homozygous: CC versus TT, OR = 1.24, 95% CI: 
1.04-1.48, P= 0.02) with moderate among-study heterogeneity 
(/=22.8%, P=0.01) (Figure 2A). Stratified analysis by etimicity 
or source of control detected no notable associations among any 
subgroup (i.e., Asian versus Caucasian populations, hospital-based 
versus population-based). We then removed one study at a time to 
explore the source of heterogeneity. It was found that exclusion of 
the study by Zienolddiny et al. [4] reduced the f value to 11.3% 
{P— 0.338), and then significance of the association no longer 
existed (Table 2). This study represented only 1 weight of the 
meta-analysis. These results suggested that this study dictates the 
among-study heterogeneity, and drives the initially observed risk 
association. While carefully reviewing Zienolddiny's study, it was 
found that this study was undertaken among smokers. Subjects in 
both case and control group were either current smokers or ex- 



smokers who have quit smoking for less than 5 years. In particular, 
mean cigarettes per day for years of smoking were 15.6 ± 8.3 and 
14.8± 6.3 for 40.4 ±12.1 and 42.3 ±7.9 years in cases and 
controls, respectively. These inclusion criteria led to a minor allele 
frequency (MAF) of 0.46 in this study that was deviated from 
MAFs of around 0.60 in Caucasian population. Therefore, the 
recruit of subjects with long smoking history might be the main 
factor accounting for the heterogeneity. 

Furthermore, four of eleven publications including six studies 
(1,151 cases and 1,084 controls) reported genotype distributions of 
rsll615 polymorphism for smokers and non-smokers. With these 
extra information based on smoking status, stratified analyses by 
smoking status showed that the association of the rsll615 
polymorphism with lung cancer risk was stronger in non-smokers 
tiian in smokers (homozygous, OR (95% CI): 2.39 (1.47-3.88)/ 
2.16 (1.41-3.30); heterozygous, 2.39 (1.44-3.95)71.77(1.21-2.60); 
dominant, 2.40 (1.49-3.85)71.94 (1.36-2.79)) without substantial 
heterogeneity (Figure 2B,Table 2). 

ERCCl rs3212986 OA polymorphism 

To investigate potential associations of ERCCl rs32 12986 
polymorphism with risks for lung cancer, 6 eligible studies with 
6,639 cases and 8,630 controls were pooled together for analysis 
[4,15,16,20,22,31]. No significant association was found (homo- 
zygous, AA versus CC, OR= 1.00, 95% CI: 0.88-1.14; hetero- 
zygous, CA versus CC, OR = 0.97, 95% CI: 0.91-1.04; dominant, 
AA7CA versus CC, OR = 0.97, 95% CI: 0.91-1.04, recessive, 
AA versus CC7CA, OR= 1.01, 95% CI: 0.89-1.15) (Figure 
3A, Table 2). In the study by Kang et al. [3 1] , genotype distribution 
in the control group did not foUow HWE. This study was included 
in the final analysis because removing it did not qualitatively alter 
the association. Moreover, no significant association was revealed 
in either the Asian or the Caucasian subgroup (Table 2). 

ERCCl rs3212961 A>C polymorphism 

To date, four eligible studies investigating the roles of ERCCl 
rs3212961 in the lung cancer risk were available [4,16,28,30]. 
These studies encompassed 1,770 cases and 1,830 controls. Pooled 
analysis failed to provide statistical evidence for a significant 
association oi ERCCl rs32 12961 polymorphism with overall lung 
cancer risk (homozygous: CC versus AA, OR = 0.87, 95% CIs: 
0.71-1.07; heterozygous: CA versus AA, OR = 0.92, 95% 
CIs:0.77-l.ll, dominant, CA/CC versus AA, OR = 0.90, 95% 
CIs: 0.76-1.07, and recessive model, CC versus AA7CA 0.91 
95% CIs: 0.79-1.06) without substantial between-study hetero- 
geneity(Table 2). Owing to the relative small number and minor 
heterogeneity of the available studies, sensitivity analysis on this 
polymorphism was not performed. 

ERCCl 3212948 G>C polymorphism 

A tagSNP ERCCl 3212948 G>C, representing common 
genetic variation in the coding region of the ERCCl gene [5], 
has increasingly received attention. Three eligible studies with 
1,537 cases and 1,835 controls were identified [5,29,30]. This 
variation was found to be significantiy associated with decreased 
lung cancer risk (heterozygous: CO versus GG, OR = 0.78, 95% 
CI: 0.67-0.90, P= 0.001; dominant: CG7CC versus GG, 
OR=0.79, 95% CI: 0.69-0.91, P = 0.001) widiout obvious 
heterogeneity (Figure 3B,Table 2). 

ERCCl rs2298881 OA polymorphism 

Four eligible studies with 4,653 cases and 6,921 controls were 
acquired for evaluating the association of another ERCCl tagSNP 
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Figure 1. Flow diagram of included studies for the associations between ERCCl polymorphisms and lung cancer risk. 

doi:10.1371/journal.pone.0097616.g001 



rs2298881 with lung cancer [20,22,28,29]. The overaU pooled 
ORs were calculated as indicated below: homozygous: AA versus 
CC OR=l.ll, 95% CI: 0.91-1.36, P=0.288; heterozygous: 
OR= 1.03, 95% CI: 0.94-1.13, P=0.53; dominant: AC/AA 
versus CC, OR= 1.04, 95% CI: 0.95-1.13), P=0.41; recessive: 
OR^l.lO, 95% CI: 0.91-1.32, P=0.322. Risk estimates 
suggested that ERCCl rs2298881 is not a risk-associated 
polymorphism in lung cancer (Figure 3C, Table 2). 

Sensitivity analyses 

As mentioned earlier, the study by Zienolddiny et al. [4] was 
removed from the final analysis for ERCCl rsll615 polymor- 
phism. Subsequendy, influence analysis showed that omitting any 
other study essentially did not alter the results, thereby confirming 
the stability of this meta-analysis. Next, the same analysis was 
performed for ERCCl rs32 12986 polymorphism, and the associ- 
ation remained unchanged, suggesting that this pooled analysis is 
also stable. 

Publication bias 

Begg's funnel plot and Egger's test were conducted to detect the 
publication bias of the meta-analysis. The shapes of the funnel plot 
for the ERCCl rsl 1615 polymorphism seemed symmetrical under 
all models (Figure 4). Nonetheless, Egger's test under homozygous 
and dominant models were significant (P= 0.03 for CC versus TT, 
P= 0.03 for TC/CC versus TT, P= 0.06 for TC vs. TT, P= 0.28 
for CC vs. TC/TT). The publication bias did not alter conclusion, 
because results in all comparison models unequivocally suggested 
that the pooled OR for the association was not statistically 
significant. Furthermore, no publication bias was detected for 
associations of rs32 12986 polymorphism with lung cancer risk 
(Figure 4), which was confirmed by Egger's test (P=0.55 for AA 
vs. CC, P= 0.16 for CA vs. CC, P= 0.08 for CA/AA vs. CC, and 
P= 0.73 for AA vs. AC/CC). 



Correlation between genotypes ERCCl polynnorphisnns 
and mRNA expression 

It was thought that ERCCl rs32 12986 and rsl 1615 may affect 
transcript stability and mRNA levels, respectively [4,5] . However, 
our results did not reveal any significant association between either 
of polymorphisms and lung cancer risk. As previously published 
[14,32], we used information from HapMap and SNPexp online 
tool to further estimate whether expression levels of ERCCl 
transcript correlate with genotypes of these two polymorphisms. 
Given subjects in our meta-analysis were from Chinese and 
Caucasian populations, we performed the correlation analyses for 
CHB (Chinese), CEU (Caucasian), and a combination of these two 
populations, respectively. Information on gen<)t\'pe of rsl 1615 
polymorphism and corresponding ERCCl mRNA expression was 
available for 32 CC (GG), 55 TC (AG), and 37 TT (AA) 
individuals. For ERCCl rs32 12986 polymorphism, there were 6 
AA, 48 AC, and 70 CC carriers with expression data (Figure 5). 
Although the expression levels oi ERCCl transcript showed a trend 
of decreasing from wide type (TT) to the homozygous variant (CC) 
genotype of rsl 1615 polymorphism in CHB population, the 
difference did not reach statistical significance. Overall, there was 
no significant difference in ERCCl transcript expression levels 
among different genotypes of rsl 1615 polymorphrsm in CHB, 
CEU, and combined populations under all the genetic models. 
Similarly, genotype oi ERCCl rs32 12986 polymorphism did not 
appear to correlate with mRNA expression level of ERCCl gene, 
either. The effect ERCCl rs32 12948 polymorphism in ERCCl 
gene expression was also explored since it was suggested to 
associate with decreased lung cancer risk. Results indicated that 
rs32 12948 was related to decreased expression levels of the ERCCl 
gene in CHB, but increased expression in CEU. (recessive model: 
CEU: P= 0.0433; CHB: P= 0.0538) (Figure SI). 
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A. ERCC1 1 9007 T>C (rs1 161 5) 
Study ID 



Caucasian 
Vogel (2004) 
Zhou (2005) 
Matullo (2006) 
Zienolddiny (2006) 
Sakoda (201 2) 
Kazma (2012) 

Subtotal (l-squared = 67.4%, p = 0.009) 

Chinese 

Yu (2003) 
Zhang (2003) 
Yin (2009) 

Deng (2011) ^ 
Yin (201 2) 

Subtotal (l-squared = 40.5%, p = 0.1 51) 

Overall (l-squared = 56.2%, p = 0.011) 
NOTE: Weights are from random effects analysis 



T 
.208 




B. 

Study ID 



SMK 

Deng (2011) 
Zhang (2008) 
Zienolddiny (2006) 

Subtotal (l-squared = 0.0%, p = 0.462) 
NSMK 

Deng (2011) 
Zhang (2008) 

Yin (200SD 

Subtotal (l-squared = 30.5%, p = 0.237) 
Overall (l-squared = 4.0%, p = 0.391) 



o 



OR (95% CI) % Weight 



1.16 (0.69,1.95) 
0.93 (0.79,1.21) 
1.18(0.68,2.07) 
2.85 (1.68,4.81) 
1.23 (0.93.1.62) 
1.02 (0.88,1.19) 
1.21 (0.97,1.51) 



0.93 (0.69,1.40) 
1.64 (0.86,3.14) 
2.17 (1.07,4.39) 
1.80 (0.88,3.70) 
0.93 (0.49,1.74) 
1.33 (0.94,1.86) 



7.39 

15.60 

6.73 

7.31 

13.55 

17.60 

68.19 



11.15 

5.48 
4.84 
4.68 
5.66 
31.81 



1.24 (1.04,1.48) 100.00 



4.81 



OR (95% CI) % Weight 



<> 



1.53 (0.47,4.99) 
1.02 (0.37,2.78) 
2.02 (1.30,3.13) 
1.77 (1.21,2.60) 



1.88(0.65,5.40) 
4.87 (1 .83,12.95) 
1.78 (0.86,3.67) 
2.39 (1.44,3.95) 

1 .98 (1 .46, 2.68) 



7.90 
12.47 
46.68 
67.06 



7.94 
6.26 
18.75 
32.94 

100.00 



.0772 
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Figure 2. Forest plots of lung cancer risk associated with the ERCC1 polymorphism. A, Forest plot of risk of lung cancer associated with the 
ERCCl rs11615 polymorphism by a homozygous model. B,Forest plot of lung cancer risk associated with the ERCC1 rs11615 polymorphism in the 
stratified analyses by smoking status. The plots of heterozygous model were shown. SMK, smoker; NSMK, non-smoker. The estimate of OR and its 
95% CI are plotted with a box and a horizontal line for each study; O represents pooled ORs and its 95% CIs. 
doi:1 0.1 371/journal.pone.009761 6.g002 



Discussion 

Accumulating studies of tlie associations between ERCCl 
polymorphisms and lung cancer risk have been conducted, but 
yielded conflicting results. The current meta-analysis examined 
whether five most commonly studied ERCCl polymorphisms 
(rsl 1615, rs3212986, rs3212961, rs3212948, and rs2298881) were 



associated with lung cancer risk. Association between ERCCl 
rsll615 or rs3212986 polymorphism and risk of lung cancer is 
seemingly biologicaUy worthy of approval, since these two 
polymorphisms have been thought to alter transcript stability 
and mRNA levels [4,5], respectively. Nevertheless, results drawn 
from the current meta-analysis did not support such an 
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A. ERCC1 8092 C>A (rs3212986) 
Study ID 



OR (95% CI) % weight 
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Figure 3. Forest plots of lung cancer risk associated with the ERCCl polymorphisms. A, Forest plot of lung cancer risk associated with the 
ERCCl rs3212986 polymorphism. The plot of dominant model was shown. B, Forest plot of lung cancer risk associated with the ERCCl rs3212948 
polymorphism. The plot of dominant model was shown. C, Forest plot of lung cancer risk associated with the ERCCl rs2298881 polymorphism. The 
plot of dominant model was shown. 
doi:1 0.1 371 /journal.pone.009761 6.g003 



association. Initially, the meta-analysis indicated a significant 
association of the ERCCl rs 11 6 1 5 polymorphism with lung cancer 
under the homozygous model with moderate corresponding 



amon g-study heterogeneity (/=22.8%, P=0.01). Nonetheless, 
influence analysis unveiled that the association mainly arose from 
one study by Zienolddiny and colleagues [4], which represented 
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A.ERCC1 19007 T>C (rs11615) 



Begg's funnel plot with pseudo 95% confidence limits 




s.e. of:logor 



B. ERCC1 8092 OA (rs3212986) 

Begg'sfunnel plot witti pseudo 95% confidence limits 
1 - 




s.e. of:logor 



Figure 4. Funnel plots of Begg's were used to detect 
publication bias on ERCC1 rs11 61 5 (recessive model) and ERCC1 
rs3212986 polymorphisms (homozygous model). No significant 
publication bias was found. Each point represents a separate study for 
the indicated association. Size of each point is proportional to its 
weight. 

doi:1 0.1 371 /journal.pone.009761 6.g004 

only 7.31% weight in the whole meta-analysis. Moreover, this 
meta-analysis indicates no statistically significant association 
between ERCCl rs32 12986 polymorphism and lung cancer risk. 

Recently, several meta-analyses have reported conflicting 
conclusion regarding to the association between the ERCCl 
rs 1 1 6 1 5 polymorphism and lung cancer risk [6-8] . Initially, a 
stratified analysis involving 4 studies of 2,279 cases and 2,808 
controls in a meta-analysis suggested no association between this 
polymorphism and lung cancer risk [7] . Cao et al. analyzed total 7 
studies [4,15,16,19,21,24,26] with 3,810 cases/4,332 controls and 
found similar results. Recently, Zhang et al. showed conflicting 
evidence, whose meta-analysis focusing on the contribution of 
ERCCl polymorphisms on overall cancer risk. Stratified analysis 
with 9 studies [4,15-17,19,21,23,24,26,27] of 4,652 cases and 
5,164 controls demonstrated that ERCCl rsll615, but not 
rs32 12986 polymorphism significandy increased lung cancer risk. 

Compared to these previous meta-analyses, the current study 
has several unique features regarding rsll615 polymorphism. 
First, three new studies [18,20,22] published in 2012 containing 



3,780 cases and 6,081 controls were included in the current 
analysis. Second, two studies [26,27] that were incorporated in 
Zhang's meta-analysis were excluded from the current analysis [6] . 
During intensive literature review, four articles conducted by Yin 
et al. [18,26-28] in the same hospital appeared to fit selection 
criterion. To eliminate the possibUity that samples might be 
repeatedly used in these studies, we ended up keeping only one 
study, which had the largest sample size, and was published most 
lately. Third, the study by Zienolddiny et al. [4] was removed from 
the final analysis due to its disproportional contribution to the risk 
estimates and among-heterogeneity. Fourth, available data were 
stratified and analyzed based on the smoking status of study 
participants. Our preliminary result suggested that the association 
between ERCCl rs 1 1 6 1 5 and lung cancer risk might be stronger in 
non-smokers than in smokers, which may be in accordance with 
the concept of genetic susceptibility that indicate individuals with 
predisposition to cancer tend to develop cancer when suffering 
low-dose hazardous exposure (e.g., second-hand smoke). However, 
due to relatively small sample size and selection bias, this result 
should be interpreted with caution. Fifth, when compared to the 
latest meta-analysis concerning ERCCl rs 1 1 6 1 5 polymorphism 
[6], samples size in this meta-analysis almost doubled. Moreover, 
In accordance with previous meta-analyses [6,8], this analysis 
confirmed that rs32 12986 polymorphism applies null effect on 
lung cancer risk. Furthermore, it was found that correlations of 
ERCCl rsll615 or rs3212986 polymorphism with corresponding 
ERCCl mRNA expression levels were negative. These results 
provide evidence supporting our findings that these two ERCCl 
polymorphisms were not associated with lung cancer susceptibility. 

ERCCl rs32 12961 polymorphism has been also well character- 
ized. With 4 studies (1,770 cases and 1,830 controls), this meta- 
analysis demonstrated no association between this polymorphism 
and lung cancer risk. We also evaluated two additional ERCCl 
tagging SNPs (rs32 12948 and rs2298881). Risk estimates suggested 
that ERCCl rs32 12948 is a protection-associated genetic variation 
in lung cancer, whereas no association was found for rs2298881. 
To our knowledge, this is the first meta-analysis that estimated 
associations of ERCCl rs32 12961, rs32 12948 and rs2298881 
polymorphisms with lung cancer risk. Nevertheless, correlation 
analysis between genotypes of rs32 12948 polymorphisms and 
mRNA expression did not agree with association study. Results for 
correlation suggested that rs32 12948 might decrease expression 
levels of the ERCCl gene in CHB, but increase the expression in 
CEU (recessive model: CEU: P= 0.0433; CHB: P= 0.0538). 
Overall, no significant correlation was found in combined CHB 
and CEU population. Several reasons may help to explain the 
discrepancy: 1) the number of available studies for rs32 12948 was 
very few, and meta-analyses with more studies are needed to 
vahdate our findings on association. 2) with only three studies, we 
were unable to perform efficient stratification analysis by ethnicity 
to investigate the association between rs32 12948 and lung cancer 
risk in either Chinese or Caucasian group. However, the results 
convey some important information: 1) ERCCl rs3212948 
polymorphism may influence mRNA expression and 2) the effect 
on mRNA expression is ethnicity-dependent. These findings may 
promote researchers to investigate biological function of ERCCl 
rs32 12948 polymorphism. 

There are some limitations in the current up-to-date meta- 
analyses. First, this meta-analysis may have been limited by the 
small number of eligible studies available. In particular, there were 
only four, three, and four publications available for rs32 12961, 
rs32 12948, rs2298882 polymorphisms, respectively, which may 
have attenuated the statistical power. Second, the missing data 
regarding cancer stage and insufficient histological details may also 
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1 . Genotyping data of ERCCI polymorphisnis and mRNA expression levels by genotypes were obtained from the HapMap 

phase ll+lll release 28 data from EBV-transfoimed lymphoblastoid eel lines ftom 270 individuals. 
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Figure 5. ERCCI mRNA expression by the genotypes of ERCCI rs11615 and rs3212986 polymorphisms. 

doi:1 0.1 371/journal.pone.009761 6.g005 



have impacts on the interpretation. Third, publication bias was 
detected by Egger's test, while analyzing the association between 
rs 1 1 6 1 ,5 polymorphism and lung cancer risk under the homozy- 
gous and dominant models. Forth, due to the deficiency of 
information on age, sex, ethnicity, and exposures (e.g., smoking, 
pack/year, and drinking), conclusions could only rely on 
unadjusted ORs, which might suffer confounding bias. Ideally, 
to perform a more precise analysis, crude ORs should be adjusted 



by potential confounders, such as age, sex, and other environ- 
mental factors after study results are pooled together. Unfortu- 
nately, in most cases, meta-analyses had to present only 
unadjusted ORs for the following reasons: 1) not all studies 
included in the meta-analysis provided adjusted ORs, 2) the ORs 
published in different studies were not adjusted by the same 
potential confounders, 3) usually only aggregate data (e.g., sex, 
age, and genotype) but not individual data are available in the 
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reports. Given these deficits, we cannot eliminate the possibility 
that the significant associations observed in the current analysis 
might result fi-om chance. Therefore, these results should be 
interpreted with caution. Further carefiiUy designed studies with 
large sample sizes of different ethnic populations are warranted to 
validate our findings. Finally, this meta-analysis has been limited 
by the failure to acquire genotype data for some planned subgroup 
analyses. Particularly, ERCCl gene polymorphism may play roles 
to different extents in smoked-related and non-smoked-related 
lung cancer; however, very few of the studies in this meta-analysis 
separately reported SNP genotype [:ounts for smokers and non- 
smokers. Future investigator should carefully list genotype counts 
for smokers and non-smokers to allow such stratified analyses to be 
done. 

Conclusion 

Collectively, with more studies incorporated, we were able to 
conduct a more precise and robust evaluation for ERCCl rsll615 
and rs32 12986 polymorphisms. We found that there was null 
association between these two polymorphisms and lung cancer 
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